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Abstract : Attempts to gain an understanding of spacecraft plasma dynamics via 

experimental investigation of the interaction between artificially 
synthesized, collisionless, flowing plasmas and laboratory test bodies date 
back to the early 1960's# In the past 25 years, a number of researchers have 
succeeded in simulating certain limited aspects of the complex spacecraft- 
space plasma interaction reasonably well* Theoretical treatments have also 
provided limited models of the phenomena. However, the available in situ data 
was fragmentary, incomplete, and unable to provide a good test for the results 
from ground based experiments and theory* Several active experiments were 
recently conducted from the space shuttle that specifically attempted to 
observe the Orbiter-ionospheric interaction* These experiments have contrib- 
uted greatly to an appreciation for the complexity of spacecraft- space plasma 
interaction but, so far, have cuiswered few questions* Therefore, even though 
the plasma dynamics of hypersonic spacecraft is fundamental to space technol- 
ogy, it remains largely an open issue* This paper provides a brief overview 
of the primary results from previous ground-based experimental investigations 
and the preliminary results of investigations conducted on the STS-3 and 
Spacelab 2 missions* In addition, several, as yet unexplained, aspects of the 
spacecraft-space plasma interaction are suggested for future research* 

1 . INTRODUCTION 

Any object placed in space will be immersed in a macroscopic ally neutral 
conglomeration of positively and negatively charged particles, as well as 
neutral particles, generally called a plasma. The space plasma is, although 
very tenuous, an important component of the space environment— of critical 
importance to many geophysical and active plasma investigations as well as to 
the environmental dynamics of large space structures such as the space 
station. However, today after three decades of space flight, there still 
remain a number of open questions and unexplained effects that require the 
attention of future research efforts. 

The physics of a body immersed in a quiescent, collisionless plasma is 
well understood: the body takes on an electric (floating) potential which 
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tends to balance the flux of charged particles to its surface so that no net 
electrical current flows. The plasma tends to shield itself from this poten- 
tial by creating a region of unequal ion and electron number density surround- 
ing the body, called the plasma sheath, in which the floating potential on the 
body is matched to the space potential of the plasma. However, when a rela- 
tive motion exists between a body and its environmental space plasma ^ an 
interaction occurs which is far more complex than the simple quiescent case. 

A redistribution of surface charge occurs on the body and the zone of disturb- 
ance in the plasma is no longer radially symmetric, regardless of body geome- 
try. When the relative motion between the body and plasma is mesosonic, as it 
is in the lower ionosphere, several characteristic processes have been found 
to occur: the plasma sheath on the frontal side of the body may be compressed 

to some extent by the directed motion of the ions; immediately behind the 
body, the more massive particles are swept out leaving a region essentially 
void of ions and neutrals; potential wells and oscillations may occur, ions 
are accelerated into the void region and ion beams and plasma waves propagate 
into, and away from, the wake. Although the electron mobility is sufficiently 
great to populate the void region fron relative velocity considerations, a 
negative space charge potential is created their presence which tends to 
impede their motion into the region. Hence, the void region that occurs in 
the wake near the body is highly depleted of all charged and neutral particles 
and forms the most intense feature of the body-plasma interaction. 

The way in which the void region is repopulated may result from a variety 
of mechanisms including the focusing of ions electric fields in the plasma 
sheath surrounding the front half of the body, ambi polar diffusion, thermal 
diffusion, the plasma expansion phenomena, and scattering by various plasma 
oscillations or instabilities. Other characteristics of hypersonic plamsa 
dynamics include the spatial extent of the interaction region, the rate at 
which the disturbances propagate outward downstream from the body, and effects 
that occur in the wake after the void region has been repopulated (Figure 
1). For example, in cases where the electrostatic focusing in the plasma 
sheath dominates, some ions may be deflected onto the trailing surfaces of the 
body and a region of significant ion number density enhancement has been 
observed to occur downstream on the wake axis at the crossing point of the 
deflected ion trajectories. 

The dominate characteristics of plasma flow interactions and the 
governing physical mechanisms depend on the various body and plasma parameters 
such as scale size, electron-to-ion temperature ratio, ion acoustic Mach 
number, and body potential — the effects of which are understood only over 
limited regions of parameter space. Moreover, the discussion, so far, has 
dealt only with simple conducting bodies and non-magnetized, collisionless 
plasmas. Clearly, other complicating factors exist, such as neutral gas 
emissions from the spacecraft, collisional effects, secondary electron emis- 
sion, solar UV, chemical reactions, multiple ion species, and magnetic fields. 
Since the relative motion between the body and the plasma is supersonic with 
respect to certain ion plasma waves, a collisionless shock wave may be expect- 
ed to occur under some circumstanes. It is further thought that secondary 
electron emission may lead to a non-monotonical matching in the plasma sheath 
of the floating potential on the body with the, generally more positive, space 
potential of the environmental plasma. Although it is recognized that these 
effects can occur in space, they have been beyond the scope of most experi- 
mental and theoretical studies, which usually treated only small scale. 
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conducting bodies of relatively simple geometry^ in collisionless / unmag- 
netized plasmas of a single ionic species (for a review see Stone, 1981a). 

Although it has long been recognized that the spacecraft- space plasma 
interaction occurs and can adversely affect spacecraft systems as well as 
scientific instrumentation (Stone et al., 1978), most space missions have 
involved single satellites which, a priori, could provide only very limited 
information. Therefore, with the possible exception of the Gemini-Agena 10 
and 11 missions, no deliberate and systematic attempt has been made to study 
the problem in space prior to the advent of the space shuttle. As a result, 
the in situ data available from the 60 *s and 70 *s are incomplete in spatial 
coverage and fragmentary in the sense that seldom were all the necessary 
measurements made (for a review see Samir, 1973). 

In the absence of definitive in situ measurements, a large number of 
theoretical and ground-based experimental studies were made during this 
period. Although much valuable information has been gained from these 
efforts, both approaches have limitatins and the field of spacecraft plasma 
dynamics, or ionospheric aerodynamics, is far from being well understood. The 
results from recent space shuttle missions have answered few questions to 
date, but rather, have contributed to a greater appreciation for the 
complexity of spacecraft-space plasma interactions. 

The intent of this paper is to touch briefly on the major contributions 
of previous experimental investigations, provide examples of the corroboration 
of the interpretation of in situ measurements by the understanding gained from 
groundbased investigations, and discuss the main features of the STS-3 and 
Spacelab 2 space shuttle investigations. 


2. MATHEMATICAL FORMULATION 

In the kinetic, formulation of plasma flow problems, since a plasma 
consists of a number of different types of particles, including electrons, 
ions, and neutral molecules, a complete statistical description of the state 
of the plasma requires a separate function, f (x,v,t), to describe the 
distribution of each constituent in six dimensional phase space at some given 
time. Each of these distribution functions is the solution of a Boltzmann 
equation of the form: 

which describes the rate of change of the distribution of the a-constituent , 

^ r in space and time. Alternatively, this equation can be viewed as a 
statement of the conservation of particles existing in an elemental volume of 
phase space, d^x d^v. 

The terra on the right-hand side of the equation. Of /9t) , accounts for 
short range, discrete collisions, which may be electrostatic ia nature 
(between two charged particles) or "hard sphere" type collisions. 
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The Lorentz force, q [i + v x S/C], results from the self-consistent 
electrostatic and magnetic fields which account for "distant collisions" of 
particles with long range forces. A collective behavior then results from 
Coulomb interactions between groups of charged particles. 

The Lorentz force is made self-consistent ty requiring it to satisfy the 
Maxwell equations: 

» X s • J || + ■l.j (2) 

7xS--|||, (3) 

and the Poisson equation: 

7*S = , (4) 

c 

where 5 is the total current flow and is the net charge density. 

The current and net charge density, in turn, depend on the distribution 
functions of the plasma constituents (solutions of the Boltzmann equations) 
through the relations 

^ " ^ext ^ ^ (5) 

ct 

and 


= E f f^(x,v,t) d^v . (6) 

a 

Equations (1) through (4), subject to the definitions (5) and (6), form 
the governing equations for plasma flow interactions. This set of partial 
differential equations is coupled and nonlinear. Therefore / a number of 
simplifying assumptions and approximations are generally made in kinetic 
treatments to obtain a tractable problem. In effect^ most experimental models 
have made many of the same simplifications. 

First/ it is generally assumed that the flow interaction exists in a 
steady state. While this obviously has a great impact on the complexity of 
the equations, its physical justification is questionable. There is no 
description of time dependent effects. There exist, however, experimental 
evidence that suggests the presence of wave particle interactions in the wake 
region. 

A second widely used assumption is that the magnetic field can be 
omitted, which reduces the Lorentz force to q i. There is some experimental 
justification for this assumption under certain conditions. This assumption 
coupled with the previous assumption that 9/^t = 0 eliminates Maxwell 
equations (2) and (3), leaving only the Poisson equation (4) that the electric 
field must satisfy. From equation (5), we see that the only currents possible 
are those resulting fron external forces, which are generally omitted. 
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The third major assumption is that the plasma is "collisionless.” This 
means thatr although long range/ collective interactions will occur and must 
be considered/ the short range/ discrete collisions occur so infrequently as 
to have a negligible effect. Since the discrete collisions are negligible/ we 
set (9f^/8t) = 0 and equation (1) becomes homogeneous (the Vlasov equation). 

Moreover/ since neutral particles interact through discrete collisions (which 
are assumed to be negligible) they can be completely neglected. We/ 
therefore/ only require time independent collisions Boltzmann • (or Vlasov) 
equations for ions and electrons; i.e../ 
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A fifth general assumption is that the flow is mesotherraal; i.e./ 

r where c ^ is the meaui thermal speed of ions or electrons and v^ 
is the orbital speedJ Since the mean thermal velocity of the ions is 
negligibly small compared to the relative motion between the body and the 
plasma/ the ions can be assumed to behave as a monoenergetic stream (no 
thermal motion). The electrons/ on the other hand/ have a mean thermal speed 
much greater than the orbital speed and maintain a Maxwellian distribution in 
the presence of a repulsive body potential. 

The boundary condition generally assumed at the body states that the body 
is an equipotential surface and that all incident charged particles are neu- 
tralized; i.e./ no charged particles are reflected from the surface. The 
potential is assumed to go to zero (plasma potential) infinitely far fron the 
body. 


3. SCALING LAWS 


The dimensionless parameters that must be invariant in order to obtain 
strict similitude between two flow interactions of different scale sizes can 
be derived formally from the governing equations (1-4). We also apply here 
the general assumptions discussed above. Making the variable substitutions; 


X = R X 
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t = T/w 


V = V u 
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f. = n F 
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7 = V/R / 
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m, = m M 
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n. = n N. , 
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the governing equations take the dimensionless form; 
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and 


4Tren R 2 
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( 9 ) 


The equations (7-9) will remain invariant if we require P = kT /e and the 
following parameter groups to remain constant: 


Z = number of charges per ion 
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Hence ^ the dimensinoless parameters, Debye ratio, R^, Ion acoustic Mach 
number, S, and normalized electric potential, arise naturally from the 
governing equations* The necessity of scaling all three parajneter groups was 
shown experimentally by Skvortsov and Nosachev (1968a) in that measurements 
taken for constant values of R^, S, and obtained at different values of 
T^, show little variation while a significant variation appears in some cases 
when only the ratios R^ and (|)^/s 2 were preserved* 

In principle, it should be possible to obtain any arbitrary combination 
of the dimensionless scaling parauneters R^, S, and by appriopriate choices 
of the physical variables n^, T^, and R^* Unfortunately, this is not the 

case since, in practice, several of these variables are subject to experiment- 
al limitations* The test body size must be smaller than the plasma stream, 
therefore, making the plasma source radius an upper bound for the body radius, 
R 9 * Further, due to the nature of ion accelerators, it is difficult to obtain 
high number densities in low energy streams* 

The above limits impose no constraint on ^ (--cl), /T ) since, in the 

laboratory, is independent of all other variables anS can be adjusted to 
any desired value by an external voltage source* Similarly, R^ (-^n R 2 /t ) 
can be made arbitrarily small and S (-V 2 /t ) arbitrarily large by maSing^R 
small and large, respectively. However, ^R^ cannot be made arbitrarily ^ 
large while making S small. Since R^ must be less than the beam radius, any 
further increase of R^ must be accomplished by increasing n^ and/or decreasing 
T^. However, this is inconsistent with small S which requires small (and 
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therefore small n^) and/or large T^« The conditions of large and small S 
are therefore mutually exclusive in the laboratory and can only be approached 
within certain practical limits. 

As a result of the practical constraints on and S, it will be possible 
to correctly scale very few of the wide range of conditions possible for 
satellites or diagnostic instruments in the ionospheric plasma 
according to the strict Vlasov scaling laws developed above. In recent years / 
however, a concept known as qualitative scaling has evolved which allows a 
considerable relaxation of the rigid Vlasov laws (Falthammar, 1974). Under 
qualitative scaling, parameters much greater (or smaller) than unity are 
required to remain so but are not required to maintain the same order of mag- 
nitude. Only paraiTieters which are of order unity must be scaled closely; i.e.. 
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Qualitative scaling has greatly extended the applicability of groundbased 
experiments to natural in situ phenomena. Further, it allows additional 
aspects of the problem, such as magnetic and temporal effects, to be included. 


4. AN OVERVIEW OF GROUNDBASED LABORATORY INVESTIGATIONS 

Most laboratory research has centered around the near to mid— wake regions 
of small bodies (1^^ ~ 1 to 10) in nonmagnetized plasma streams. The primary 
results for this case will be discussed below. For a discussion of the 
upstream and far-wake disturbance, the reader is referred to the work of 
Fournier (1971), Hester and Sonin ( 1970a, b), and Woodroffe and Sonin (1974). 
The dynamics of magnetized plasmas has been studied by Astrelin et al. (1973) 
and Bogasheohento et al. (1971). 

• Simulation of the Spacecraft-Ionospheric Interaction 
( 1 ) The Disturbance Envelope 

The envelope of the zone of disturbance, defined by the boundary 
between freestream conditions (Ji/J^o = 1) and disturbed flow (J^/J^o ^1), 
depends on two factors; the initial width of the disturbance at^the largest 
cross section of the test body, and the rate at which the disturbance propa- 
gates away from the Z-axis as it moves downstream. The initial radial extent 
of the disturbance, defined by the sum of the test body radius and the sheath 
thickness, was found to increase in proportion to the ion acoustic Mach number 

and the negative body potential as (Figure 2). It can also be 

expected to increase with the Debye length. The propagation of the disturb- 
ance boundary away from the wake axis was found to define a Mach cone 
(Skrortsov and Nosachev, 1968b, and Stone et al., 1978) based on the ion 
acoustic Mach number, S. This result is in agreement with several theoretical 
treatments, including those of Rand (1960a,b) and Maslennikov and Sigov (1965, 
1967, 1969) vdiich predict a Mach cone structure for bodies with a small 
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potential ^ . The rarefaction wave, which is the most spatially extensive 

characteristic of the disturbance, was found to decrease the ambient ion 

current density by as much as a factor of three at distances as great as 

2(S*R ) downstream (Stone et al., 1978). 
o 

The above conclusions are based on the parameter range = 4 to 
6, S ^ 11, and - -3.8 to -47. Measurements by Hester and Sonin ( 1969a, b; 
1970) show the existence of pseudowaves (streams of ions deflected across the 
wake axis by the sheath fields) for >>S. These ion streams may overrun the 
rarefaction wave, under certain conditions, and extend the zone of disturbance 
beyond the mach cone. 

(2) Ion Trajectory Focusing by the Plasma Sheath 

The focusing of ion streams onto the wake axis by the electric 
field existing in the plasma sheath surrounding a test body was inferred in 
early studies by the presence of diverging wave-like structures in the far- 
wake region (Hester and Sonin, 1970a, b and Stone et al., 1972, 1974). More 
recent, direct vector measurements of deflected ion streams show their angle 
of attack to be proportional to in the near wake for small scale bodies 
(Stone, 1981b). ^ 


The vector ion flow measurements also revealed a "bunching” of 
the ion trajectories at the radial boundaries of the ion void region, which 
can be seen in the current density profiles of Figure 3. This effect was not 
discussed in early theoretical or experimental studies (Stone, 1981a), 
although it is apparent in the theoretical results calculated by Maslennikov 
and Sigov (1969). 

The present experimental data are not sufficient to reveal the 
physical mechanism which produces the observed ion trajectory "bunching” in 
the near wake. Neither do the calculations of Maslennikov and Sigov allow an 
explanation, and the the effect does not even occur in other theoretical 
treatments such as the one by Call (1969), which predicts a "fanning out” of 
the deflected ion trajectories. The near-wake ion trajectory grouping may be 
produced either by a collective effect on the ions (possibly a result of 
instabilities set up by the large density gradient at the void boundary) or by 
a nonmonotonic potential gradient in the plasma sheath. 

(3) The Axial Ion Peak 

Early investigations by Hall, Kemp, and Sellen (1964) clearly 
show the ion void in the near-wake region and an axial ion peak for a 
spherical test body. The ion current density was measured at a number of 
stations across the wake for a wide range of body potentials, revealing a 
distinct dependence of the axial ion peak on More detailed measurements 

of the ion peak, including both transverse and axial profiles for a variety of 
potentials, were published a year later by Clayden and Hurdle (1966). These 
measurements, in addition to showing a dependence on show the axial ion 

peak to rise rapidly behind the body and trail off slowly, extending more than 
20 R^ downstream. Similar results were found for a sphere and a conical body 
oriented with its apex into the flow. 
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Converging ion streams at the boundaries of the wake void region 
were found to create the initial ion peak on the wake axis ( Stone ^ 1981c). 

The position of this peak was found to depend on S, R , and as shown in 
Figure 4. This result agrees closely in its S and dependences with the 
theoretical predictions by Martin (1974), A second peak may also created on 
the wake axis further downstream. The two types of axial ion peaks are norm- 
ally superpositioned for small and cannot be distinguished from each other. 
However, they tend to separate at highly negative values, indicating that a 
second casual mechanisms with a different dep)endence on $, may be involved— 
possibly the collisionless plasma expansion phenomenon. “ 

The height and width of the cucial ion peak at the location of 
its maximum amplitude were also investigated (Stone, 1981c). The maximum peak 

height for spherical test bodies was found to be proportional to [S/|$ I] as 

b* 

shown in Figure 5. The peak width (normalized the test body radius, R^) 

was found to depend only on | | ( pj^g^^e 6). This dependence is taken to 

represent a balance between the momentum the particles obtain due to deflec- 
toward the Z— axis, produced by (p.r and the magnitude of the space charge 
potential barrier on the Z-axis, which is proportional to kT . (Note 
that = e(j.^ATg.) ® 

It was also found that the nature of the axial ion peak depends 
strongly on the cross-sectional geometry of the test body. This can be 
explained simply by the behavior of the plasma sheath, which is directly pro 

portional to 9^S). For small and large R^, the sheath is thin compared 

to the test body dimensions and conforms closely to its geometry. For a test 

body having a square cross section, the ion streams were found to be deflected 
onto lines that are orthogonal at the Z— axis and approximately a body diameter 
in length. This produces a wider peak of lower amplitude. As becomes 
large or R^ becomes small, the sheath becomes relatively thick, expands away 
from the body and acquires a more spherical shape. Hence, the ions are 

more toward a point on the Z-axis and the axial ion peak structure 
will take on a behavior more characteristic of spherical test bodies. 

The effects of test body geometry were studied in a preliminary 
manner and indicate that the wake of a geometrically complex body may be 
explained in terms of a linear superposition of the wakes of its different, 
simple, geometric constituents (Stone, 1981a). 

(4) Deflection of Ion Trajectories in the Mid-Wake Region 

The deflection of ion trajectories in the mid-wake region was 
first predicted theoretially by Maslennikov and Sigov (1969) and later by Call 
(1969) although not to the same degree due to the limitations imposed by his 
flux tube technique. In both cases, the ion trajectories were found to be 
deflected away from the Z-axis by the positive space charge potential 
associated wih the mid-wake axial ion peak, this effect was inferred by 
Hester and Sonin (1970) from the nature of the diverging wave— like structure 
they observed in the far wake. More recent vector ion flow measurements show 
that ion streams exist within the plasma wake with angles of inclination to 
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the Z-axis smaller than the geometric angle defined by the radial extent of 
the test body (Stone, 1981a)* Since ion streams obviously cannot pass through 
the test body and do not originate on its rear surface, this is taken as clear 
evidence that the streams, initially deflected toward the Z-axis, subsequently 
underwent an additional deflection away from the Z-axis somewhere downstream — 
presumably at the location of the axial ion peak* 

(5) Effects of Ion Thermal Motion 

An extensive experimental study of this effect was carried out 
by Fournier and Pigache (1975) and it has been studied theoretically by Taylor 
(1967), Gurevich, et al* (1969), and Fournier (1971)* While it is clear that 
the general effect of ion thermal motion is to diffuse the detailed wake 
structure discussed above, the more quantitative question as to how effective 
this diffusion process is, or at what value the structure vanishes, has 

not been satisfactorily resolved. The calculations by Fournier for a long 
cylindrical body show that a nonmonotonic n^ distribution continued to exist 
in the wake for < -2*75 when - 2 and for < -6 when = 1* 

Moreover, the experimental observations by Fournier and Pigache (1975) show 
the peak structure to completely vanish only for small and < 1 

(Figure 7)* It appears that the opposing effects of and are such 

that for small the axial ion peak vanishes for - 1/ but can be 

recreated by a sufficiently negative value of $^^* We may, therefore, conclude 
that in the ionosphere, where for - 2 and ^be axial ion peak 

can be expected to occur to some extent* This conclusion is supported by the 
clear presence of an axial electron peak (which in laboratory studies is 
smaller than the ion peak) in the wakes of the Ariel I satellite and its 
spherical ion probe (Henderson and Samir, 1967). 

(6) Effects of Large Values 

Laboratory and theoretical studies have been limited to 
values less than 50, which do not approach the R^ range of large space 
platforms. However, several parametric trends have been established that may 
be cautiously extrapolated to describe large bodies in Earth orbit* 

It was pointed out above that the amplitude of the axial ion 
peak significantly depends only on S and and, further, that the peak width 
(normalized by body radius) depends only on <^^* If these observations, made 
over a relatively small range of R^, can be extrapolated to large R^ values, 
then, on this basis, the axial ion peak would be expected to maintain an 
approximately constant width (relative to the body radius) and amplitude as R^ 
increased arbitrarily. 

This conclusion is incomplete, however, without considering the 
effects of ion thermal motion* The tendency of random motion to spread out 
and diminish the wake structure can be expected to increase with the distance 
traveled by the ions, and hence with R^* Therefore, it becomes doubtful that 
the detailed wake structure discussed above would be observed for very large 
bodies in the ionosphere at floating potential* However, if the body poten- 
tial became elevated (such as may occur in the case of the Space Shuttle 
orbiter when charged particle accelerators are fired or as the result of 
(v X B) • £ potentials on very large structures such as the Space Station) the 
wake structure may appear as a result of the opposing effect of $ shown by 
Fournier* ^ 
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(8) Magnetic Field Effect 

The effects of the geomagnetic field on spacecraft plasma 
dynamics have been included in only a few laboratory investigations; e.g., 
Bogashchento et al. (1971) and Astrelin et al. (1973). In particular, it was 
found that the net result of a parallel magnetic field is the generation of 
standing axial ion current density oscillations along the wake axis with a 
period proportional to (Za)^^/2rV ). Evaluating these results for typical 
ionospheric conditions shows that such oscillations would be very small in 
cunplitude and that the period would extend fetr beyond the mid-wake zone. 
Therefore, it is concluded that the omission of the geomagnetic field in 
studies of the near— and mid— wake regions of small to rnedixm sized bodies is 
justifiable. 

2. Process Simulation Experiments 

As stated in section 3, process scaling involves the philosophy of 
* 3 'i*J-itative scaling to study an individual physical process that may be one of 
many involved in a complete phenomenon. Activity in this area has picked up 
considerably in recent years. We will briefly discuss three examples: the 

investigation of collisionless plasma expansion across a strong density 
9^^®<3ient, the possible coupling between widely separated current sources in 
the magnetized ionospheric plasma, and electron heating of the near wake. 

(1) Expansion of a Collisionless Plasma; 

The collisionless expansion of a large reservoir of plasma 
across a strong density gradient is greatly Influenced ty an electric field 
created at the expansion front by high speed electrons separating from the 
massive ion component. The ions are, in turn, accelerated to high velocities 
(several times the ion acoustic speed) by the field, which is maintained by a 
continual replenishment of fast electrons at the front from the plasma reser- 
voir . A recent review of the phenomenon and its possible occurrence in space 
plasmas is given ^ Samir et al. (1983). A laboratory process simulation 
investigation of collisionless plasma expansion is in progress at MSFC (Wright 
et al., 1985; 1986, and several other institutions (Chan et al., 1984). 

Figure 8 is a schematic of the basic process as depicted by the 
self-consistent theory. Note that ion acceleration is constant over the 
period for which the self-consistent equations remain valid. An experimental 
study by Wright et al. (1985) was conducted in a steady-state plasma flow with 
the density gradient created by a plate oriented perpendicular to the flow. 
Figure 9 provides a comparison of experimental data with the expansion front 
velocity predicted by the theory. It is apparent that with proper scaling of 
S and this process must be simular to the process(s) responsible for 
filling in the near-wake void of orbiting spacecraft; particularly in the case 
of large scale spacecraft. 

(2) A Tethered Satellite-Electron Beam Current System 

Stenzel and Urratia (1986) have investigated the current flow 
between a field-aligned electron beam and an electrode collecting return 
current on a different flux tube in a large laboratory magnetoplasraa. This 
investigation has revealed the existence of anomalous cross-field currents 
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that shunt the field aligned current system, and temporal current disruptions. 
If these processes reflect a valid picture of large scale current systems in 
the earth's ionosphere, this study will be of vital importance to investiga- 
tions planned with long, conducting tethered satellites (e.g., the NASA TSS-1 
mission) or charged particle beam experiments on large scale structures such 
as the Space Station. 

(3) Enhancement of Electron Temperature in the Near-Wake 

An apparent electron temperature enhancement, which coincided 
spatially with the ion void region, was reported by Samir et al., 1974 (Figure 
10). The heating of electrons may be partially explained by the effects 
expected to occur in the presence of a potential well, or they may result from 
a two-stream instability produced by an interaction between fast and slow 
moving ions. 


5. COMPARISON OF LABORATORY AND IN SITU RESULTS 

A number of laboratory experiments have allowed direct quantitative 
comparison with, and sometime a better understanding of, in situ 
measurements. Here we consider of few examples. 

Several observations of an elevated electron temperature in the wake of 
ionospheric satellites have been reported (Samir and Wrenn, 1972; and Troy et 
al., 1975). This effect was investigated in the laboratory by Illiano and 
Storey (1974) and by Samir et al., (1974), the latter revealing electron 
temperature enhancements in the wake of a test body in a collisionless 
streaming plasma as discussed above in Section 4.2.3. The enhancement ranged 
up to 200 percent of the ambient stream's temperature and was confined to the 
void region of the plasma wake - in complete agreement witlwith the in siti: 
observations. The Samir et al. study did not establish the mechanism for 
electron heating in plasma wakes, but suggested it may result from a potential 
well in the void region or from wave particle interactions (instabilities) 
within the strong density gradient at the wake boundary. 

Stone and Samir (1981) report on a comparison between laboratory wake 
experiments in which ion focusing by plasma sheath electric fields were 
studied (Stone, 1981b) , and in situ observations of structure in the wakes of 
the Ariel 1 satellite and its spherical ion probe (Henderson and Samir, 

1967). The laboratory results show the ion current peaks, observed at only a 
one axial distance in the in situ wake, to be part of an extended complex 
structure such as shown in Figure 3 (Stone and Samir, 1981). Moreover, the 
effects of body potential on the wake structure observed and quantified in the 
laboratory (discussed in Section 4. 1.2-3) explain the similarity observed 
between the wakes of the satellite and the small, but negatively biased, 
spherical probe. 

A second example presented by Stone and Samir (1981) shows that the 
ram/ wake current ratio data, obtained from the AE-C satellite (as a function 
of plasma composition, electron temperature and satellite potential) can be 
collapsed to a single curve using the body-potential and ion acoustic Mach 
number dependences established in the laboratory (section 4.1.3). This 
example also shows that dimensionless parcimeters must be calculated using 
specific ionic mass and concentration values rather than average values. 


138 



6. INITIAL RESULTS FROM SPACE SHUTTLE MISSIONS (STS-3 AND SL-2) 

By far/ the most elaborate insitu investigations of spacecraft plasma 
electrodynamics to date have been conducted from the space shuttle. It was 
anticipated that the more detailed shuttle-borne experiments would complement 
existing in situ data and provide a sufficient data base for resolving several 
spacecraft— space plasma interaction issues* However the orbiter environment 
proved far more complex than anticipated and the main contribution to date has 
been a greater appreciation for the complexity of spacecraft- space plasma 
interactions* 

The third space shuttle mission (STS— 3) provided the first opportunity to 
measure the Orbiters plasma and field environment* This was accomplished on 
mission days three and four by manuevering the Plasma Diagnostic Package (POP) 
up to 15 m above the Orbiter payload bay with the Remote Manipulator System 
(RMS) as/ for example/ shown in Figure 11* Differential vector measurements 
of ion flow direction/ current density/ and energy were made during this 
period with a Differential Ion Flux Probe (DIFP) (Stone et al*, 1985). These 
measurements revealed the existence of secondary ion streams in the vicinity 
of the Orbiter at high angles of attack as great as 50® with respect to the 
ram direction and from 10 to 40% of the ram current intensity (see Figure 
12). The source or generating mechanism for these high inclination secondary 
ion streams was not identified* However/ their energy was close to that of 
the ram ions and it was concluded that they were not of geophysical origin 
(Stone et al./ 1983). 

The existence of ion streams in the disturbed plasma surrounding an 
orbiting body is not surprising; in fact/ the existence of ion streams was 
inferred by Henderson and Samir (1967) and has been studied in the laboratory/ 
as discussed in section 4*1*2* In all previous caseS/ the ion streams were 
associated with the wake region downstream fr<xn the satellite or test body and 
such streams were anticipated in the wake of the Orbiter. However/ the 
secondary streams observed near the Orbiter during the STS-3 mission were 
totally unexpected in that they were measured when the POP was not in the 
Orbiter* s wake and/ in some cases/ when it was extended upstream from the 
Orbiter* 

Futher analysis revealed several additional effects. Not only were high 
inclination ion streams typically observed/ but the ion ram current direction 
did not, in general, correspond closely to the direction of the orbital 
velocity vector (see Figure 12) . Moreover, in one case in which the PDP was 
extended above the Orbiter and oriented such that the DIFP faced directly into 
one of the secondary streams, the stream vanished (Stone et al*, 1986). The 
ion current density of both the Ram and secondary streams was also found to be 
directly proportional to the neutral particle density as shown in Figure 13 
(Stone et al./ 1986). 

These observations should be considered in the context of wave 
measurements that revealed the existence of broad band (30 to 178 kH^) 
electrostatic noise in the vicinity of the Orbiter (Shawhan et al*/ 1984) and 
the observation of higher than normal ambient ionospheric plasma densities 
(Raitt et al./ 1984; Siskind et al., 1984). 
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The above observations, taken together, are evidence that the Orbiter 
travels within a neutral gas cloud that results fron outgasing, thruster 
burns, waste dumps, etc., and that the presence of this gas cloud significant- 
ly affects the way in which the Orbiter interacts with the ambient mangeto- 
plasma. The proportionality of the ram and secondary ion stream intensities 
to the density neutral particles very near the Orbiter suggest a very 
effective ionization mechanism. The vanishing of the secondary streams at 
about 10 m ahead of the Orbiter suggests that the interaction with the 
ionosphere may be confined within an envelope that extends on the order of 10 m 
in the forward direction but presumably trails out to greater distances down 
stream. If an electric field exists at this gas cloud- ionosphere interface, 
the ambient ions would be deflected such that the ramed ion current would be 
skewed from the Orbital velocity vector at all points except in the Orbiter* s 
XY-plane, where the envelope would be normal, and hence the electric field 
parallel, to the velocity vector. This is in agreement with the trend of the 
data as shown in Figures 11 and 12. The motion of the secondary ion streams 
through the background plasma represents a large source of free energy that 
may be expected to generate instabilities. Recently Hwang et al. (1986) have 
developed a model which shows that ion streams in a background ionospheric 
plasma can generate electrostatic noise over the oberved spectrum. 

7. CONCLUSIONS AND OPEN ISSUES 

The application of laboratory plasma physics to space plasma physics and 
technology has proven successful in addresssing a number of issues over the 
past 25 years; particularly in addressing, through qualitative simulation, the 
overall interaction of small ionospheric satellites and instriiments with the 
ionospheric plasma (for a review see Stone, 1981a). In the past few years, 
process simulation has been used to address several specific aspects of 
natural physical phenomenon, such as plasma expansion (Wright et al., 1986) 
and the bahavior of current systems in magnetosplasmas (Stenzel and Urratia, 
1986). The conclusions reached in these studes have not been verified by 
in situ data yet and their usefulness remains to be determined - although the 
preliminary analysis of Space lab 2 data suggests the existence of a process 
very much like the collisionless plasma expansion process studied in 
laboratory plasmas by Wright et al. 

Through the combined efforts of laboratory and in situ investigations 
together with theoretical treatments, the physics of geometrically simple 
conducting bodies of 1 to 50 Debye lengths scale size is reasonably well 
understood — although certain details remain to be determined; e.g. the effects 
of a magnetic field on the far wake and the mechanisms producing changes in 
the electron distribution. This is not true, however, either for geometric- 
ally complex bodies (although previous studies (Stone, 1981a) suggest a linear 
superpositioning of the wakes of the simple geometric constituents, this must 
be more universally established), bodies with complex surface characteristics, 
bodies with an associated neutral gas cloud, large scale size bodies (in 
excess of 1,000 Debye lengths), or bodies of any type with an elevated elec- 
tric potential. Future investigations of spacecraft plasma electrodynamics 
should therefore center around five issues; i.e., 

(1) The effects of body scale size, particularly for large scale 

bodies (with respect to both the Debye length and cyclotron radius). 
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(2) The physics of a neutral gas cloud in a hypersonic magnetoplasnia« 

(3) The effect of body potentials well in excess of the ionization 

potentials for neutral constituents. 

(4) The effects of complex spacecraft geometry. 

(5) The effects of complex spacecraft surface characteristics. 

Moreover, these issues are of such complexity that their resolution will 
require laboratory investigations (that offer the advantage of comparatively 
low co.st and fast turnaround) and theoretical treatments (that make use of 
recent gains in computer technology) , as well as deliberate and systematic 
in situ experiments designed to provide data capable of corroborating the 
ground-based results. The previous in situ results from small ionospheric 
satellites were not systematic and too limited in scope, while the experiments 
conducted on the space shuttle involved an extremely complex body. Clearly, 
the required in situ experiments must be systematic, closely controlled, and use 
cleaner test bodies (in terms of geometric complexity, surface 
characteristics, gas emissions, and EMI). An understanding of the fundamental 
physics of hypersonic spcecraft plasma dynamics will be important to the 
P^^ctical application and operation of large space platforms and should not be 
taken lightly in planning wide usage of the space station. 
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Figure 1. Ion behavior within 
disturbed zone* After Stone 
[1981a] . 
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Figure 2. Dependence of the sheath 
width at the critical cross-section 
of a test body {x/\ ) on body poten- 
tial and ion acoustic Mach number. 
Open circles represent data for a 
sphere and open triangles for a 
parallel cylinder. After Stone et 
al. 


Figure 3. Ion current density 
profiles downstream fron a 
conducting sphere (R^ = 3 cm) for 
= 0.8^ S = M, and = -5, data 
after Stone et al. [1972]. 
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Figure 3. Schematic of plasma Figure 9. Ion velocity at the 

expansion into a vacuum. (a) expansion front vs time. O, labor- 

initial condition, ( b) Evolution of atory measurements, — , theoretical 

density, (c) Evolution of ion veloc- model. After Wright et al. [1986] . 

ity according to the self-similar 
treatment. After Samir et al. 

[1983] . 
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Figure 10. Transverse profiles of 
[Tg (wake)/T^^j downstream from a 
conducting sphere for = 1200 K® , 
count n = 7.5 x 10^/cm^, - 5.3 

eV. After Samir et al. [1974]. 



Figure 11. POP track on Julian day 
85 for the period 16:48:40 to 
16:51:05. DIFP normal is indicated 
by n. Dashed lines indicate the 
inferred boundary of the interaction 
region. After Stone et al. [1986]. 
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Figure 12. Angles- of- attack of the 
ram ion current (open circles) and 
secondary ion steam (closed circles) 
for the period indicated in Fig* 11* 
Solid line is the angle between the 
DIFP normal and the orbital velocity 
vector. After Stone et al. [1986] • 


Figure 13. Ion current density vs 
pressure (neutral particle density) 
for the period indicated in Fig. 

11* After Stone et al. [1986]. 
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